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ABSTRACT: We aim at reducing the uncertain- 
ties inherent in the analysis of the topological 
structure by using scale controlled smoothing and 
observables independent on the "microscopic" de- 
scription of the instanton ensemble. 

Investigating the role of instantons in the QCD 
vacuum by lattice simulations is a challenging 
problem. Instantons are involved in the J7a(1) 
symmetry breaking (cf the Witten - Veneziano 
formula which relates quenched topological sus- 
ceptibility and 77' mass) and in chiral symmetry 
breaking (via zero modes of the Dirac operator), 
and they lead to dynamical effects at intermedi- 
ate distances. Global properties like susceptibil- 
ity and charge distributions afford direct tests, 
while the features of the local structure (ensem- 
ble, size and distance distributions) are ingredi- 
ents for dynamical models To recall, an R'* 
(anti-)instanton (A)I of size p located at x = 
has charge (q) and action (s) densities: 
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Q = y d^xq{x), J '^''2; s(a;) = IQI 5o 

s{x) — So \q{x)\, So — Stt^, Q = iinteger. 

"Superpositions" of N Fs or A's lead to higher 
minima of the action: S — NSq. However pairs 
I- A are not minima and, e. g., under unrestricted 
cooling they decay, the sooner the smaller their 



action 2 5o - < 2 So is, where 

depends in particular on the "overlap" u = {pi + 
PA)/diA, diA being the I - A distance. 

In analyzing the topology by lattice methods we 
need to deal with the roughness of the Monte 
Carlo configurations at short scales and to iden- 
tify the physically relevant topological structure. 
UV lattice artifacts (dislocations) and close I-A 
pairs, indistinguishable from short range density 
fluctuations, have small action and can be easily 
produced. To smooth out in a controlled way high 
frequency fluctuations in Monte Carlo configura- 
tions the method of Restricted Improved Cooling 
(RIC) has been developed in |^ as a gauge invari- 
ant low pass filter. RIC introduces a parameter 
directly related to a physical scale above which 
fluctuations will be preserved. Since it uses an 
"improved" action RIC has as fixed point non- 
trivial classical configurations. But since it pre- 
serves all structures above the chosen scale it also 
retains, e. g., I-A pairs with overlap smaller than 
a threshold depending on the smoothing scale. 

The identification of the uncovered topologi- 
cal structure poses special problems. Usually a 
"microscopic" description in terms of I's and A's 
is attempted. This, however, becomes increas- 
ingly ill defined at small scales, particularly if 
close pairs abound. Therefore it is interesting to 
avoid the necessity of such a description and ob- 
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tain the phenomenological parameters appearing 
in the instanton models from observables which 
can be directly measured on the lattice. 

We here suggest to study properties of instan- 
tons in lattice simulations using observables ob- 
jectively defined in continuum QCD. Rather than 
investigating instantons microscopically by in- 
specting lattice field configurations, we shall con- 
sider "macroscopic" observables which in a pure 
instanton vacuum can directly be related to prop- 
erties of the instanton ensemble. One such class of 
quantities are ratios of VEV's of chirally odd op- 
erators, which are purely non-perturbative quan- 
tities, and (at least in a dilute instanton vacuum) 
are independent on the instanton density and can 
directly be related to moments of the instanton 
size distribution. Assuming the latter to be a well 
defined property of the topological structure, we 
expect these ratios to remain approximately con- 
stant under cooling, at least after achieving a cer- 
tain, minimal degree of smoothing. Using RIG we 
can rephrase this question in terms of the physical 
scale of the relevant topological fluctuations. 

Scale Controlled Smoothing. Cooling ||] is an 
iterative, local minimization of the action, which 
proceeds sweep-wise and can be defined to con- 
verge onto (non-trivial) classical configurations. 
It acts as a diffusion process with the length scale 
of smoothing growing like the square root of the 
number of iterations 

Restricted Improved Cooling (RIC) is a 
scale controlled smoothing procedure involving 
two ingredients: 

a) - Improved minimization action with practi- 
cally scale invariant instanton solutions, and 

b) - Restriction of cooling to allow a certain 
amount of (Euclidean) "energy" above the mini- 
mum, homogeneously distributed over the lattice. 

Our improved action 0, H is correct to order 



O(a^) and is completely flat for instanton sizes 
larger than po ^ 2.3a, below which it drops - see 
Fig. ^ It then follows from a) that: 

- Below the "dislocation threshold" po short range 
topological structure is eliminated. Note that 
Pa Q in continuum. 

- Above Po, instantons are stable to cooling. 

- The correspondingly improved charge density 
leads to a charge approaching an integer already 
after a few cooling sweeps and stable thereafter. 
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Figure 1. Improved (5Li) and Wilson action vs 
instanton size given in lattice units - these are 
SU(2) results, but they are representative also for 
SU(3). 

The restriction b) is introduced as the following 
modification of the local updating rule for a link: 

U iff A = TiiWW^ - {UW^)^) > 6^ (2) 

where W is the staple connected to U and V the 
group projection of W (for SU(2) F = VF/||Ty||). 
It naturally leads to gradual saturation and even- 
tual stop of cooling after a number of sweeps de- 
pending on 5. RIC thus acts as a frequency filter 
for the field fluctuations. Since A is the square 
of the lattice equations of motion it has a contin- 
uum limit and therefore the cooling parameter 5 
can be related to a physical scale. 
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no cooling 


144(2) 





.5754 


.01 


411.3 


.27(5) 


176(3) 


8 


.9846 


.04 


290.8 


.34(6) 


178(3) 


16 


.9888 


.04 


145.4 


.41(7) 


181(3) 


30 


.9939 


.05 


51.41 


.67(9) 


183(3) 


67 


.9976 


.07 



Table 1 

RIC results for SU(3). 
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Figure 2. Mass gap M{t) — a Lvs lattice distance 
t {L: spatial lattice size). The horizontal band 
indicates standard results for this lattice. 



For this we calculate the string tension a from 
correlations of spatial Polyakov loops at separa- 
tion t in time, measured on configurations RI- 
cooled with 5. We determine rc{5) as the the 
minimal distance t at which the string tension is 
recovered, taking the first cooling curve as refer- 
ence. See Fig. |^; for details and for SU(2) results 
see 1^. The data presented here concern SU(3) 
and have been obtained on a 12'', pbc lattice at 
/3 = 5.85 (o = 0.135fm), using 900 configurations 
separated by 600 sweeps (after 50000 thermaliza- 
tion sweeps). The line in Fig. || corresponds to 
re((5) « 2.3 (5-1/3. 

The restriction b) does not affect significantly 




Figure 3. RIC smoothing scale rc{5) vs 5 ^Z'^. 



Figure 4. Polyakov loop distribution in complex 
plane. Left: no cooling, right: RIC {5 = 290.8). 

the behavior of single instantons under RIC but 
does affect the behavior of I-A pairs. There is a 
well defined relation between A and the I-A "in- 
teraction" S\^^ = 167r2 - S*^^ g. As a result 
RIC with given 5 preserves pairs depending on 
their interaction, hence on their overlap. Gener- 
ally, physics on scales larger than Tc is expected to 
remain unaffected by RIC. How RIC affects small 
and large scales is also illustrated by the values 
of the Plaquette and Polyakov loop (|P|) in Ta- 
ble 1. The distribution of the latter after cooling 
stays compatible with confinement (see Fig. ||). 
An extended topology analysis for SU(2) using 
RIC has been provided in here we add some 
further results for SU(3). 

The topological charge distribution is Gaussian 
and is stable under coohng - see Fig. ||. In 
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Figure 5. Topological charge distribution from 
RIC at (5 = 411.3 (squares), 290.8 (diamonds, left 
shift) and 6 — 145.4 (crosses, right shift). The 
tendency to depletion in the |Q| = 1 sector with 
increasing smoothing is due to the periodic b.c. 




Figure 6. Density correlations vs lattice distance 
(binned) for no cooling and for RIC with the 
first three values for S in Table 1 (curves from 
up to down). The right plot shows the region of 
the minimum, enlarged, for the same cooled data 
(curves from left to right). 



Table 1 we give the topological susceptibility of 
SU(3) for various i5 (the phenomenological expec- 
tation is x^/'' 180 MeV). On Fig. | we plot the 
SU(3) density correlations {q{x) q{0)) at various 
values of S. Due to the improved charge density 
even non-cooled data can be obtained, they show 
however strong UV renormalization effects. No- 
tice that {q{x)q{0)) < for disjoint supports |^]. 
With progressing smoothing this effect disappears 
indicating gradual loss of support properties. 

From the RIC analysis of SU(2) 1^ it appeared 
that one can speak of a typical size, but that 
a "microscopic" description of the instanton en- 
semble is problematic since topological excita- 
tions cannot be separated at small scales from 
short range fluctuations. We observed at small rc 
a strongly growing 1(A) density, with more and 
more peaks of alternating charge and increasingly 
large overlap showing up at distances in the Vc 
range. This may also affect the size determina- 
tion and explain the diferences in the detail of 
the size distribution found in the literature ^J. 
To avoid the uncertainties introduced by this sit- 
uation in the determination of the instanton size 



we look for an independent determination of the 
typical size using "macroscopic" observables. 

"Macroscopic" observables for topology. 

The relevance of instantons to chiral symmetry 
breaking is due to the fact that a single / {A) 
induces a localized (i.e., normalizable) zero mode 
of the Dirac operator of definite chirality: 

zV(A,(^))$±(a:) =0, 75 $±(a;) = ± $±(x). (3) 

Instanton models usually proceed from a sim- 
ple picture of the Yang-Mills vacuum as a dilute 
"medium" of /'s and A's: 

p^n-^l, or p/i?<l (4) 

where p is the average size, R the average distance 
and n — N/V the density of the instantons. The 
ratio p/R can be used as a small parameter to 
classify non-perturbative effects generated by the 
medium of instantons In this picture 

chiral symmetry breaking can be understood as 
a collective effect involving all instantons in the 
ensemble [|ll|, |l2j (this can be seen as a con- 
sequence of the cumulative effect of the chirally 
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odd 't Hooft vertices of individual instantons, re- 
sulting in the appearance of a dynamical quark 
mass) . The chiral order parameter is then Jl^ : 

Generally, the VEV's of any chirally odd op- 
erator can serve as order parameter for chiral 
symmetry breaking. Since they are purely non- 
perturbative quantities which acquire a non-zero 
value only because of the spontaneous breaking of 
chiral symmetry they provide good probes of the 
instanton effects. So, for instance, in an instanton 
medium we have: 

(Vi T[F] T^)^-p-'' {-p^nf^ ^-p-" (p/R) ' (6) 

where TlF] is a function of the gauge fields, T = 
1; 75i CTp^ a chirally odd Dirac matrix and d is the 
mass dimension of the operator. Note that while 
these VEV's depend differently on the instanton 
size, they all show the same dependence on the 
instanton density as the usual quark condensate. 
Therefore in ratios of such VEV's the dependence 
on the instanton density cancels: 

X={i,T[F]Tij)/{{pij)^p-''^' (7) 

For example one finds in the large-iVc limit, 
where all instantons are of size p — p, [ p^ : 

{^F^,a^,ib)/{i^^Jj)=Ap-^ (8) 

In the instanton vacuum such ratios can thus di- 
rectly be related to properties of the instanton 
size distribution, with no reference to the density. 
(A non-zero density is needed only for numera- 
tor and denominator to be non-zero individually.) 
Therefore we expect quantities of the type (|^) to 
be approximately "invariant" under cooling. 

We plan to measure the cooling behavior of ratios 
of the type ^ using RIG. Since we can control 
the smoothing scale we can test the dependence of 
our derivation on the assumption of diluteness by 



observing the saturation of the ratios, and then 
obtain an estimation of the typical instanton size 
independently on a microscopic description of the 
ensemble. (A lattice simulation for (H) has been 
performed with limited statistics in [p^.) 
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